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A STUDY OF INTERNAL AND DISTRIBUTED DAMPING 
FOR VIBRATING TURBOMACHINER BLADES 


The principal purpose of the research was to study internal 
and distributed damping as possible methods for reducing the 
vibration response of turbomach inary blades, and to develop 
theoretical methods for analyzing the damped vibration. Efforts 
were focussed on continuous, rather than localized, damping. 

Thus, the widely used damping mechanisms such as rubbing between 
shrouds, platforms, or in blade roots, which already have received 
wide-spread attention, were not considered. 

This work, which was terminated after the first third of what 
was to be a more comprehensive study, was limited to the forced 
vibration analysis of blade-like models subjected to viscous and 
material (i.e., hysteretic) damping. Both distributed and 
concentrated exciting forces were considered. It was necessary to 
make theoretical extensions of the well-known Ritz-Galerkin 
analytical methods to deal with the out-of-phase response of 
damped systems. The developed procedure was successfully applied 
to cantilever beam, plate and shallow shell representation of 
blades, and is particularly useful for making parametric studies 
of the effects of damping upon blade vibration response. 

The theoretical development and application is summarized in 
two publications. One is a paper which was presented at the 
Vibration Damping Workshop in Long Beach, California, February 27- 
29, 1984 ("Extensions of the Ritz-Galerkin Method for the Forced, 
Damped Vibrations of Structural Elements," Flight Dynamics 



Laboratory Report AFWAL-TR-84-3064, pp. EE-1 to EE-22). Another 
is the Ph.D. dissertation of Mr. T. H. Young (Ohio State 
University, Department of Engineering Mechanics) which is 
currently being completed. 

A literature search was made to uncover previously published 
work dealing with the forced vibration responses of continuous 
systems with damping. Approximately 150 relevant references were 
found. A review article summarizing the majority of these 
references was begun, but was dropped when project funding was 
di scontinued. 



EXTENSIONS OF THE RITZ-GALERKIN METHOD FOR THE FORCED, 
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ABSTRACT 

The classical method for analyzing the forced vibrations of structural 
elements such as beams, plates and shells is to express the displacements as 
superpositions of the responses of the free vibration modes. This is only 
possible for those relatively few problems where exact eigenfunction solutions 
exist, and often only with considerable difficulty. Ritz-Galerkin methods are 
widely used to obtain approximate solutions for free undamped, vibration prob- 
lems. The present paper demonstrates how these same methods may be used 
straightforwardly to analyze forced vibrations with damping. This is done 
directly without requiring the free vibration eigenfunctions. 

The Galerkin method fias been shown previously to be an effective tech- 
nique for these types of problems. The Ritz method has the further advantage 
of not needing to satisfy the force-type boundary conditions, which is par- 
ticularly important for plates and shells. But proper functionals represent- 
ing the forcing and damping terms must be developed, and this is done in the 
present work. 

The present paper considers two types of damping— viscous and material 
(hystereticl . Both distributed and concentrated exciting forces are treated. 
Numerical results are obtained for cantilevered beams and rectangular plates. 
Studies showing the rates of convergence- of the solutions are made. In the 
case of the cantilever beam, approximate solutions from the present methods 
are compared with the exact solutions. 
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I . INTRODUCTION 


The Rayleigh-Ritz methods are probably the best known and most widely 
used methods for analyzing the free, undamped vibrations of more complicated 
structural elements, such as beams Ce.g., variable thickness), plates and 
shells. The Rayleigh method, which has been in existence for over a century 

[1] , takes advantage of the fact that a system vibrating in one of its normal 
modes, interchanges its energy completely between potential and kinetic forms. 
That is, the maximum potential energy (V max ) equals the maximum kinetic energy 
CTmay) which occurs in a cycle of vibration. Assuming a trial function for 
the mode shape, one which satisfies at least the geometric boundary conditions 
of the structural element, and assuming simple harmonic motion, setting T^x 
equal to V max yields the vibration frequency. Unless one is fortunate enough 
to have assumed the exact eigenfunction (i.e., mode shape) of free vibration 
for the trial function, the resulting frequency is too high, and is conse- 
quently an upper bound on the exact solution. 

In 1908 Ritz 12] improved upon the Rayleigh procedure by assuming a set 
of admissible trial functions, each having independent amplitude coefficients. 
He was able to show that a closest upper bound for the frequency could be 
achieved by minimizing the functional Tmax" v max with respect to the coeffi- 
cients. Indeed, he applied this method to the completely free square plate 

[2] for which no exact solution is possible. Since then literally hundreds 
of references may be found which apply the method to free vibration problems 
(cf.[3,4]]. In some cases such as shells, an exact solution may exist, but 
it is so complicated to use that a Ritz procedure may be employed more easily 
to obtain accurate results (e.g. , [5-9] ) . But in most cases no exact solution 
exists, and an approximate, properly convergent technique such as the Ritz 
method becomes essential./ 

The classical method for solving forced vibration problems for structural 
elements is to expand both the forcing function and the displacement response 
in terms of the free vibration eigenfunctions . This requires first solving 
the free vibration problem and finding an orthogonal set of eigenfunctions. 
Further, the forced vibration response typically requires the tedious (and 
computationally expensive) evaluation of integrals of products of the eigen- 
functions with themselves and the loading functions. If the eigenfunctions 
are complicated (for example, combinations of the regular and modified Bessel- 
functions in the case of a circular plate) the integrals may have to be evalu- 
ated numerically. In many cases (e.g., beams) large roundoff errors are easi- 
ly generated during the procedure. 

It can be shown that the method of Galerkin [10] is equivalent to the 
Ritz method, and is some ways more general (it may be applied even when no 
energy functionals exist for the problem) . This equivalent method has also 
been frequently used to analyze the free vibrations of structural elements 
such as beams, plates and shells (cf.,[3,4]). Several years ago a method of 
applying the Galerkin method to forced vibration problems of continuous 
systems having viscous damping was demonstrated [11,12] . The method is a 
direct one; that is, the forced response is found without first having to 
solve for the free vibration eigenfunctions. The method was demonstrated for 
the one-dimensional, second order differential equation of the vibrating 



string, and the two-dimensional, fourth order differential equation of .the 
vibrating plate. 

However, the Galerkin method has one disadvantage in comparison with the 
Ritz method. Both methods require satisfaction of only the geometric boundary 
conditions. However, if the force-type boundary conditions are not also satis- 
fied, the Galerkin method requires the use of additional terms. For a plate 
having free edges, for example, these additional terms consist of two line 
integrals taken along each free edge which must be added to the orthogonalizing 
integrals. In the case of a shell, four line integrals would have to be added 
for each free edge. Thus for a cantilevered plate or shell of rectangular 
planform, the Galerkin method becomes quite cumbersome. 

The Ritz method has been demonstrated in recent years to be an efficient, 
accurate technique for the analysis of free vibrations of cantilevered shells 
having rectangular planform 113-16] . Such configurations may be used for tur- 
bomachinery blade vibration studies . The method is particularly useful for 
preliminary design stages, or for making parameter studies, where finite ele- 
ment methods have been found to be very costly and unreliable. 

The present paper summarizes recent research in extending the Ritz method 
to forced vibration problems where damping is present. This eliminates the 
aforementioned disadvantages of the Galerkin method when free edges are en- 
countered. The crux of the problem is to find a dissipation functional ac- 
counting for the damping forces which are 90 degrees out of phase with the 
exciting, restoring and inertia forces of the system, which may be added to 
the kinetic, strain and load potential energies in a suitable manner to give 
the correct solution. The extension is demonstrated in the present work for 
systems having either viscous or material damping, distributed or concentrated 
exciting forces, and one-dimensional (beam) or two-dimensional (plate) struc- 
tural elements. The convergence of the extended method to exact solutions is 
demonstrated for a set of cantilever beam problems. A similar method was de- 
veloped by Siu and Bert £17] for laminated composite plates having material 
damping, subjected to distributed forcing pressure. 


II. BEAM ANALYSIS - VISCOUS DAMPING 


2.1 Distributed Forcing Function 

Consider first a beam of length l subjected to a distributed transverse 
load q (force/unit length) which varies sinusoidally with time, 

q(x,t) =q(x)e int (2.1) 


where (2 is the forcing frequency and x is the coordinate measured along the 
length. The kinetic energy during vibratory motion is 


_ 1 r , 3w % 1 
- 2 / 


dx 


( 2 . 2 ) 


where w is the transverse displacement, p is the mass per unit length and A is 
the cross-sectional area. The strain energy of the beam due to bending defor- 
mation is 



(2.3) 



EK^d* 

3x 


where E is the modulus of elasticity and I is the second moment of the cross- 
sectional area with respect to the neutral axis of bending. 

Define a dissipation functional V by 

V » j I c f| wdx (2.4) 

o 

where c is the viscous damping coefficient. It should be noted that V differs 
from the well-known dissipation function of Rayleigh [1]. The work done by 
the exciting force- is 

£ 

W = / qwdx (2.5) 

o 


To apply the Ritz method, 
may be expressed as 

w(x,t) = W(x)e^ t 

- f C <j> (x)e int 
j=l J J 


assume that the vibratory displacement w(x,t) 

( 2 . 6 ) 

(2.7) 


where the 4>-j (x) are trial* functions which satisfy at least the geometric boun- 
dary conditions (zero displacement and/or zero slope) which exist for the 'beam, 
and the Cj are complex coefficients which may be separated into real and imagi- 
nary parts as 

C. = C* - iC 1 (2.8) 

3 3 3 

R 

where i=/-l . The real part (C.) may be regarded as the vector component of 
the response in phase with the^exciting force, and the imaginary part (Cj) is 
the response component which lags the exciting force by 90 degrees. 


Define further a generalization of the functional T mav -V mav used in free, 
undamped vibration analysis by the Ritz method. That is, let 


L = (T -V ) - (V -W ) (2.9) 

max max max max max 


where the separate terms on the right-hand- side of Eq. (2.9) are the maximum 
values of the functionals previously given by Eqs. (2.2) through (2.5) , ob- 
tained by substituting Eq. (2.6) into them and setting |e2iQt|-i. That is, 


T 

max 


2 £ 

fl / ,2. 

— J pAW dx 


(2.10a) 



(2.10b) 


V 

max 



V 

max 



W 

max 


£ 

/ q W dx 
o 


(2.10c) 
(2 .10d) 


The Ritz method is applied by substituting Eq. (.7) into Eqs.(2.10) and by 
using the minimizing equations 

3L 

“ 557 s - * o ( j =1 <r 2 , . . . J) (2.11) 

which results in a set of J linear simultaneous equations in the unknown coef- 
ficients Cj . Setting the real and imaginary parts of each equation equal to 
zero yields a set of 2J equations in the unknowns and C* , with coupling 
between the Cj and the C? coefficients arising from the damping. The right- 
hand sides of the equations evolve from the forcing function (q) . Solution of 
this set of equations completely determines the forced, damped response. 

2.2 Example . Cantilever Beam with Uniform Pressure 

Consider the uniform^ homogeneous, cantilever beam of length £ clamped 
at the end x=0 and free at the end x= £ (see Fig. 1) , subjected to a uniformly 
distributed pressure 

q(x,t) = q o e li2t (2.12) 



Figure 1. Cantilever beam with uniformly distributed forcing pressure . 



where q Q is a constant. Let the trial functions to be used in Eq. (2.7) be 
expressed as simple algebraic polynomials, J 

<j>..(x) = x j (j=2,3,...,J) (2.13) 

which, by beginning with j=2, guarantees that the boundary conditions at the 
clamped edge 

3w 

w(0,t) = r-(0,t) = 0 (2.14) 


are satisfied exactly. Substituting Eqs.(2.12) and (2.13) into Eqs.(2.7), 

(2.9) and (2.10), carrying out the minimization indicated by Eqs.(2.11) and 
the necessary integrations over x, and separating the real and imaginary parts, 
yields the following set of equations: 


TeI j ( j-l) k(k- 

1) _ pAQ 2 1 _R 

r cQ 1 r 1 

L 4 j+k_3 

j+JM-lJ C j + 

[j+k+l J 3 

- r cn l C R + 

[ei j (j-l)k(k-l) 

PM 2 1. 

[ j+k+l J c j 

Li 4 j+k “ 3 

j+k+l J 

These 2 (J-l) 

equations may be 

solved for 


o 

k+1 


,J) 


(2.15a) 

(2.15b) 
„ I 


comparison, and yielded Eqs.(2.15) identically. 

The response of the beam at a typical point may now be obtained by adding 
the in-phase and out-of-phase components separately, and combining them vec- 
torially. That is. 


w (x, t) = W(x)e l(fit ' <f) 

with the amplitude of the response given by 


(2.16) 


W(x) = 


r r J 

r R , 

2 

V I, 

2 

l C -cj> (x) 
_j=2 3 J _ 

+ 

l C.^.(x) 

L j=2 3 3 J 



and the phase angle lag by 

J 


(J) (x) = tan 


l cV(x) 

-1 j=2 3 J 

r R rf> 


(in radians) 


(2.17) 


(2.18) 



2.3 Exact Solution of the Previous Problem 


The preceding problem of the cantilever beam subjected to uniformly dis- 
tributed, sinusoidally oscillating exciting pressure has an exact solution. 
The equation of motion for the problem is 


trT 3 4 w . , 3 2 w 3w , , iQt 
EI — + pA — + c = q (x) e 

3x 4 3t 


(2.19) 


Assume a solution to Eq.(2.19) in the form 


w(x,t) = l X m (x)T m (t) (2.20) 

m=l 

where X m (x) is a typical eigenfunction of the free, undamped vibration problem. 
That is , 

X (x) = cosha E-cosaE-y(sinhaE-sinaE) (2.21) 

m m m m m m 


where £ = x/£ and cMl is the nondimensional frequency parameter given by 


(a £) 2 = ai i 2 J 
m m V EI 

which are the eigenvalues of the equation 

cosa £ • cosha £ = -I 
m m 


( 2 . 22 ) 


(2.23) 


The coefficient y is determined from 
m 

cosha £+cosa £ 

■ m m 

^m sina £+sinha £ 
m m 


(2.24) 


Values of a £ and y are given in Table 1. 
m m 

Substituting Eq.(2.20) into Eq.(2.19), multiplying through by a typical 
eigenfunction x > integrating over the length, using the orthogonality rela- 
tionships 

a a 4 v 

/ V. = 0 = ! ~TT V 12 - 

o o dx 

(when m^n) and the equation of motion for free vibrations 


d \ 2 

EI r - pAw X„ = 0 

4 m m 

dx 


(m=l,2,. . .) 


(2.26) 


Table 1 . Eigenfunction parameters for clamped- free beams . 


m 

°m £ 

Ym 

1 

1.8751041 

0.7340955 

2 

4.6940911 

1.01846644 

3 

7.8547574 

0.99922450 

4 

10.9955407 

1.00003355 

5 

14.1371684 

0.99999355 

• • • 
CO 

(2m-l)ir/2 

• • • 
1 


yields the equation for the response of the mth mode: 

pAT + cT + pAu) 2 T = q_e^‘^ t 
m m mm m 

where the dots represent time derivatives, and where 
l 

/ q o*m dx 




/ ^m 2 ** 


( 111 = 1 , 2 . . .) 


The solution to Eq.(2.27) is 

T (t) = (A -iB )e lQt 
m mm 


(m=l , 2 . . .) 


where A and B are given in nondimensional form as 
m m 

6 .2 

A = -2.(1 - JL ) 
m A (1 w 2 ' 
m 

.5 „ 2 

B = -£- S-) 

m A c til 
cm m 


4 = 

a) cm m 

m 


(2.27) 


(2.28) 


(2.29) 


(2.30) 


Here 6 is the static amplitude of the mth mode (i.e., when fl/cu =0) and c 

m . m m cm 

is the critical damping constant of the mth mode. In detail. 


6 = 
m 


Sn 


pAtu 


- I c = 2pAw 
2 cm m 


m 


(2.31) 


The total response at a typical point on the beam is given by 
w(x,t) = Ce 1 


(2.32) 


where C = 


oo 


00 

£ Vm (x> 

2 + 

l B X m (x) 
. mm 

|_ m=l J 


,m=l 


(2.33) 


and , 

<p = tan — 


l Vm (x) 

m=l 


(in radians) 


(2.34) 


It is observed that C and <f> are both functions of x. 

The classical solution procedure shown above is exact in that explicit 
expressions for responses of the free vibration modes (i.e, Eqs.(2.30)) may be 
written, even though the total response requires infinite summations of terms. 
The crux of the problem ie the evaluation of the integrals of Eq. (2.28) . For 
this purpose the tables calculated by Young and Felgar [18] are useful. 

2.4 Numerical Results for the Uniform Pressure Loading 

Table 2 presents numerical results for the amplitude response C/<5 Q of a 
cantilever beam subjected to uniform pressure varying sinusoidally with time, 
where 6 Q is the static displacement (0=0) of the point considered. In this 
case the viscous damping is small (c/c cl =0.01) , and response at the free end 
(x=£) is evaluated. The frequency ratio (0/oj^) is varied from 0 to 2, and 
data are given also for excitation at the second and third natural frequencies 
(Q/w^ = 6.267 and 17.547). Solutions using the Ritz method are obtained using 
2, 4 x and 7 polynomial terms of the type given by Eq. (2.13) , yielding 4, 8 and 
14 simultaneous equations from Eqs.(2.15). Comparison is made with the 
results of the exact solution, described in Section 2.3. 


It was found that, to at least five significant figure accuracy, the 
exact solution is the same as the 7-term Ritz solution for all except at 

the third resonance and that (again, except near the third resonance) a 4- 
term solution would be sufficiently accurate for engineering purposes. If 9 
terms are used, five significant figure agreement with the amplitude at the 
third resonant point is also found. Results for the response at the middle 
of the beam (x=0.5£) yielded nearly identical comparisons. 



Table 2. Amplitude response C/6 0 at the end (x=£) of a cantilever beam; uni- 
formly distributed load, viscous damping (c/c c j_=0 .01) . 


Number of trial functions in Ritz method 
fl Exact 


oil 2 4 7 Solution 


0.0 

0.5 

0.9 

0.99 

1.0 

1.01 

1.1 

1.5 

2.0 

6.267 

17.547 


1.00000 

1.33134 

5.06759 

28.6556 

45.8860 

44.5984 

5.04832 

0.82658 

0.34731 

0.03837 

0.00001 


1.00000 

1.33759 

5.29616 

36.0870 

50.6697 

35.5744 

4.81325 

0.82473 

0.35276 

1.27647 

0.02056 


1.00000 

1.33759 

5.29622 

36.0896 

50.6694 

35.5715 

4.81315 

0.82472 

0.35276 

4.48085 

0.83682* 


same 

It 

II 

ll 

It 

If 

II 

n 

n 

n 

0.93829 


*9-term solution yields 0.93829 


The static solution (fl/to-^O) in this case is exactly expressible as a 
polynomial of 4th degree. Thus a 3-term Ritz solution would be exact at this 
frequency, while the "exact" method must laboriously sum up an infinite set of 
eigenfunctions. Table 2 shows the 2-term solution to be also exact at the 
particular point x=£. At x=0.5il, it would yield 0.941176, instead of unity. 

Interestingly enough, the exact solution procedure gave numerical diffi- 
culty in obtaining accurately converged results for Table 2, even with the 
use of double-precision (i.e., 16 significant figure) arithmetic. This was 
largely for large a l, the Y m are very close to unity, and cosh a m 5 and 
sinha m £ are nearly identical. However, the Ritz method proceeded straight- 
forwardly without difficulties. 

2.5 Point Loading 


Consider next the case of excitation by a concentrated transverse force 
F^ = F-.e^^ t acting at an arbitrary point x=x^ along the beam. Solution of the 
problem by the Ritz method is simple and straightforward. Instead of Eq. (2.5) 
the work done by the exciting force is 


W = F 1 w(x 1 ) 


(2.35) 


and the appropriate functional for use with the Ritz method is 


0/ ‘ = 
max 


FjWfx^ 


(2.36) 


instead of Eq. (2.10d) . When two or more point loads are present, or if a point 
load is acting in addition to a distributed load, then superposition applies. 






In the case of a point load acting at the free end of a cantilever beam, the 
right-hand- side of Eq. (2.15a) is simply replaced by F . 

To obtain an "exact" solution to the problem in terms of free vibration 
inodes when a concentrated force is present, either of two methods may be used. 
One method would represent the concentrated force as- a Dirac-delta function 
for q(x,t) , and expand this function in terms of the eigenfunctions as in Eq. 
(2.28). However, this procedure would be extremely tedious, for it would re- 
quire many terms in the series (Eq.(2.20)) to determine a reasonably accurate 
representation. 


Another exact procedure introduces a change of variables. For example, 
consider the end loaded beam. The equation of motion is homogeneous: 


4 2 

_ ,_3 w. .3 w ^ 3w 
El — — + pA — — .+ Ct— = 0 

3x 4 3t 2 


(2.37) 


But one of the boundary conditions is nonhomogeneous; i.e., 

w(0, t) = r— '(0 , t) = t) = 0 

3x 3x 2 


(2.38a) 


EI-— r(£ ,t) = F e lJit 
3x J 


(2.38b) 


Let w (x , t ) be replaced by 

w(x,t) = v(x,t) + g(x)e 


mt 


(2.39) 


where g(x) is chosen so that 


g ( 0 ) = |3-(0) = 2_£(jD = o 
3X 3x 2 


(2.40a) 


ei L^.( £ ) = F x 


(2.40b) 


3x' 


Then g(x) is found to be 

, . F 1 2 + F 1 3 

* (x) = - iiI X + 6 iT X 


(2.41) 


Substituting Eqs.(2.39) and (2.41) into (2.37) transforms the problem into one 
having a nonhomogeneous equation of motion, with homogeneous boundary condi- 
tions ; viz . 


3 4 v 

El— + pA— 
3x 4 3 1 


= (pAft 2 -iCft) g (x) e ifit 


+ 


(2.42a) 



(2.42b) 


v( 0 ,t) = fj( 0 ,t) = ^U,t) 

3x 


.3 

3 v 

= o 
3x J 


and the solution proceeds as in Section 2.3. 

III. BEAM ANALYSIS - MATERIAL DAMPING 


3.1 Incorporation of. Complex Stiffness 

In the case of material damping (also called "structural" or "hysteretic" 
damping) , it is possible to use an equivalent viscous damping representation 
under certain conditions, the viscous damping coefficient being chosen so as to 
dissipate an equivalent amount of energy per vibratory cycle. The most impor- 
tant condition to be met is that, for the given forcing function and the range 
of frequency ratio under consideration, a single normal mode is strongly domi- 
nant among the modes excited. While it is possible to represent any single 
mode response adequately by equivalent viscous damping, especially in the vi- 
cinity of resonance and if the damping is reasonably small, the individual 
modes respond differently. Therefore, strong normal mode coupling makes the 
representation less accurate. 

The preferred method of treating material damping is by means of a complex 
modulus of elasticity (cf. [19]). That is, let the modulus be 

E* = E(1 + in) (3.1) 

where i=/d and n is the loss factor. 

* 

To apply the Ritz method to a forced vibration problem with material damp- 
ing, the energy dissipation due to friction is combined with the elastic strain 
energy in complex form, and the functional to be minimized is 

L = T - V* + W (3.2) 

max max max max 


In the case of a beam. 


for example. 


V* is given by 
max 


V* 

max 


1 3 
- i f E(l+in) dx 


3x 


(3.3) 


whereas Tg^^ and Cl ) mav remain as given previously by Eqs. (2.10a) and (2.10d) , 
respectively . 

The exact solution uses the equation of motion expressed in complex form, 
which is 

4 2 

.. .,3 w t ,3w , , int 

E(l+in)I — 7 + PA — — = q(x) e 

3x 4 3t 


(3.4) 



Assuming that the displacement (w) and the transverse loading (q) can both be 
represented as summations of the eigenfunctions of free vibration, a procedure 
similar to that employed in Section 2.3 may be followed, leading to the solu- 
tion forms given by Eqs.(2.29) and (2.32) , except that in the case of material 

damping the coefficients A and B become 

m m 


A 

m 



B 

m 





A 

s 


(1 - 


ft 


0 ) 


2 2 


+ n 


m 


(3.5) 


with 6^ given by Eq. (2.31) , and the amplitude and phase angle again given by 
Eqs. (2.33) and (2.34), respectively. 

3.2 Example. Cantilever Beam with Uniform Pressure 


The uniform, homogeneous, cantilever beam subjected to sinusoidally vary- 
ing, uniform pressure, as described previously in Section 2.2, is again taken 
as an example; however, in this case the damping is hysteretic. Using alge- 
braic polynomial trial functions as in Eq. (2.13) with the Ritz method, the 
following set of equations arises for the determination of the coefficients Cj 
and C? associated with the real and imaginary parts of the solution: 

J t 

y f [" El j(j-l)k(k-l) pAfl 2 1 R I" El j (j-l)k(k-l) "I Ii _ q o 

ji 2 { [7 j^*=3 j^rj c j + [ n 7 jTk=3 J c j } - k^r (3 - 6a 

j ( j-1) k(k-l) 1 R [El j ( j-l)k(k-l) _ pAft 2 
j+k-3 I ] I 4 j+k-3 j+k+1 

Application of the Galerkin method to the problem yielded Eqs. (3.6) identi- 
cally. The displacement w(x,t) of an arbitrary point is again given by Eqs. 
(2.16), (2.17) and (2.18). 

Numerical results for the nondimensional amplitude ratio C/S Q at the free 
end of the beam as a function of the frequency ratio are presented in 

Table 3 for a loss factor representative of many metals (n=0.001) , and in 
Table 4 for a large loss factor (q=0.1). Once again it is seen that the Ritz 
method converges rapidly to the exact solution as polynomial terms are added 
to the displacements. A 7-term solution is accurate to 6 significant figure 
accuracy except in the vicinity of the third resonance. 

In studying Tables 3 and 4, the data is seen to exhibit certain charac- 
teristics similar to those for a single degree of free system having material 
damping (e.g., (20], p. 230). That is, 


C^} = 0 


(3.6b) 




Table 3. Amplitude response C/5^ at the end (x=!) of a 
cantilever beam; uniformly distributed load, material 
damping (n=0.001). 


fl - 

Trial functions- 

-Ritz method 

Exact 

“1 

4 

7 

Solution 

0.0 

0.999999 

0.999999 

same 

0.5 

1.33770 

1.33770 

It 

0.9 

5.31980 

5.31986 

H 

0.99 

50.8387 

50.8463 

H 

1.0 

1013.39 

1013.39 

n 

1.01 

50.3766 

50.3686 

If 

1.1 

4.83954 

4.83944 

If 

1.5 

0.834967 

0.824959 

If 

2.0 

0.352793 

0.352792 

ti 

6.267 

1.32170 

14.3001 

i« 

17.547 

0.025578 

0.925293* 

1.06942 

*9-term solution yields 1 

.06942 



* 


Table 4. 
cantilever 
damping ( n 

Amplitude response C/5 Q at the end 
beam; uniformly distributed load, 
=0.1) . 

(x=£.) of a 
material 

a 

Trial functions-Ritz method 

Exact 

ui 

4 

7 

Solution 

0.0 

0.995037 

0.995037 

same 

0.5 

1.32593 

1.32593 

•f 

0.9 

4.70715 

4.70719 

it 

0.99 

9.93494 

9.93498 

•i 

1.0 

10.13252 

10.13249 

n 

1.01 

9.93658 

9.93649 

n 

1.1 

4.36873 

4.36865 

it 

1.5 

0.822066 

0.822058 

it 

2.0 

0.352364 

0.352362 

n 

6.267 

0.148371 

0.145056 

ti 

17.547 

0.009427 

0.010733 

0.010730 




(1) At n/u) =0, C/6 is not unity (as in viscous damping) , but C/6 = 

1 2 1/2 ° o 

tl/(l+n )) ' . 

(2) Damping does not shift the first resonant peak - it remains at £2/u) = 
1 . 1 

(3) The first resonant amplitude ratio (C/6 ) is very nearly equal to 

l/n. ° 

That C/6 q at fl/(u^=0 is not exactly l/n is due to the slight influence of the 
higher modes. Comparing Tables 3 and 4 , it is seen that increased material 
damping causes greater contribution from the first mode at the second and 
third resonant frequencies, as might be expected. 

IV. PLATE ANALYSIS 

4.1 Ritz Method 

Except for the fortuitous case of a plate having at least its two oppo-r 
site sides simply supported, no exact solution for the free vibration eigen- 
functions is possible [3]. Nevertheless, the Ritz method may be generalized 
from the preceding beam analysis to deal with the forced vibrations of plates 
having arbitrary edge conditions, including- point supports. 

To be somewhat specific, let us consider a rectangular plate having di- 
mensions a x b, such as the cantilever shown in Fig. 2. The exciting force 


' y 



Figure 2. Cantilever plate with concentrated exciting force. 

may be either distributed or concentrated. A sinusoidally varying distributed 
load will, in general, be of the form 

i£2t 


q(x,y,t) = q(x,y) e 


(4.1) 



where now q is force per unit area, whereas a point load is of the form 


F = Fe iflt 


(4.2) 


The transverse displacement is taken as 
w(x,y,t) = 


(4.3) 


J K 


* l l C -ik < * l -ik (x ' y) 
j=l k=l 


iflt 


(4.4) 


where the 4> j satisfy the geometric boundary conditions exactly, and the Cj^ 
are complex coefficients given by 


„ R .„I 

C., = C .. - iC., 
:k 3k 3k 


(4.5) 


The functionals for the plate, which are the generalizations of the beam 
functionals given by Eqs. (2.10) , are: 


fl 2 r a r b/2 2 
T = —- / / phW dxdy 

2 o -b/2 


(4.6a) 


max 2 


a b/2 

= T / / D{(V 2 w) 2 - 2(l-v) 


o -b/2 


3 2 w a 2 w 
L3x 2 3y 2 


2 

(— ) 
l 3x3y ; 


dxdy 


(4.6b) 


• n a b/2 

V f f 

^ 2 o -b/2 


c W dxdy 


(4.6c) 


a b/2 

W = / / q W: dxdy 

max J J 


o -b/2 


(4.6d) 


In these expressions the flexural rigidity (D) is 
3 


D = 


Eh' 


12(1 


-v 2 ) 


(4.7) 


with h being the plate thickness and v being Poisson's ratio. In the case of 
material damping, the integrand for v max is simply multiplied by (1+in) . 

The functional L is constructed as in Eq. (2.9) and minimized as 
max 


3L 

max 

3C. . 
3k 


= 0 


( 3 — If • • / k 1 , • . • , K) 


(4.8) 



R I 

yielding a set of 2JK simultaneous equations in the coefficients C., and C. . 
Solution of these equations determines the forced, damped response, J 


4.2 Example. Cantilever Square Plate with Point Load 


To demonstrate the Ritz method, a cantilever plate of square planform 
(a/b=l) is chosen, and is excited by a concentrated force F acting at the 
center of the outer edge, the force varying sinusoidally with time (Fig. 2) . 
Assume viscous damping affects the plate uniformly (i.e., c is a constant). 
Let Poisson's ratio (v) be 0.3. 


A set of polynomials for use in Eq. (4.4) may be chosen as 
<f>j k (x,y) = x^y k (j=2,3...J; k=0,l...K) 

* which satisfy the geometric boundary conditions 


w(0,y,t) = — (0,y,t) = 0 


(4.9) 


(4.10) 


exactly. Moreover, for the present example problem where the exciting force 
preserves the geometric symmetry already present, the antisymmetric modes are 
not excited, and the index k may be taken as k=0,2,4.. . with no loss of gener- 
ality. 


Accurate numerical results were obtained by taking 6 terms in the x- 

direction (i.e., j=2,3,...,7) and 4 terms in the y-direction (i.e., k=0,2,4,6). 

This yields 6 x 4 = 24 polynomial terms and, hence, 24 modal degrees of 

freedom to the problem. Considering the m-phase and out-of-phase components 

of motions, the resulting problem consists of solving 2 x 24 = 48 simultaneous 

equations for the C?. and C* . 

;ik ]k 

Numerical results for a small damping ratio (c/c c ^=0.01) are shown in 
Table 5. Here, as in the beam analysis c c ^ (=2pho) 1 ) rs the critical damping 
coefficient for the first free vibration mode of the square cantilever plate. 
In Table 5 the nondimensional displacements WD/Fa 1 2 of three significant points 
are given, as functions of These points are the plate center (x=0.5a, 

y=0) , the point of loading (x=a,y=0) and the corner points (x=a,y=+0.5b) . In 
addition to the first resonant frequency (f2/u)^=l) , data is also given for the 
next three frequencies of symmetric modes (0/ui-|_=6.1317, 7.8307 and 15.6165) . 
(i)3/u)l=7.8307 corresponds to a chordwise bending free vibration mode, whereas 
the other three are spanwise (i.e., flapwise) bending modes (cf . , [3,13]). 

Comparing Table 5 for the plate with Table 2 for the beam, both having 
the same damping ratio (c/c c ^=0.01) , one observes that the amplitude responses 
at the second and third spanwise bending frequencies is considerably greater 
for the plate than the beam. This is mainly due to two differences: 

(1) The concentrated end force excites the higher spanwise bending modes 

more than a uniformly distributed force. 
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Table 5. Displacements WD/Fa of a cantilever plate with a 

sinusoidally varying point load (P) at the center of the outer 

edge; viscous damping (c/c =0.01). 

Cl 


a 

“l 

Plate center 
?=°- s '£ = ° 

Load point 

£= i, £= o 

a b 

Corner point 

0.0 

0.111 

0.360 

0.329 

0.5 

.150 

.473 

.440 

0.9 

.601 

1.798 

1.744 

0.99 

4.153 

12.209 

11.982 

1.0 

5.755 

16.887 

16.597 

1.01 

3.992 

11.690 

11.507 

1.1 

.548 

1.574 

1.572 

1.5 

.096 

.247 

.268 

2.0 

.043 

.089 

.114 

6.1317 

1.714 

3.997 

1.843 

7.8307 

1.193 

1.254 

2.865 

15.6165 

.736 

1.258 

1.307 


(2) The concentrated force is not along the entire free end of the plate, 
but only at its center, thereby causing a significant contribution 
of the first chordwise bending mode. 

It is interesting to Observe the change in chordwise bending effects in 
Table 5 as n/u^ is swept from zero to the fourth natural frequency. At ft/u^ = 

0 one has static loading, and the plate exhibits significant anticlastic curva- 
ture; that is, the center of the free edge deflects more than the comer 
points. The anticlastic curvature is greater here than for a uniformly loaded 
square cantilever [22], as might be expected. As is increased, the curva- 

ture remains anticlastic until slightly past the first resonance (Q/u^ = 1.1). 
The curvature is strongly anticlastic at the second resonance, but reverses 
itself for the third and fourth resonances. 

Table 6 gives the results for the same problem when the damping coeffi- . 
cient is ten times greater (c/c cl =0.1). The responses at the resonant peaks 
are all very similar to those with small damping, except with amplitudes es- 
sentially one-tenth as large, except for the fourth resonance. 

V. CONCLUDING REMARKS 

An extension of the Ritz method has been developed for analyzing the 
forced vibrations of structural elements such as beams, plates and shells in 
a straightforward and efficient manner. Both viscous and material damping 
may be accounted for. The exciting forces may be distributed and/or concen- 
trated. Convergence to the exact solutions for a set of cantilever beam prob 
lems was demonstrated. 
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Table 6. Displacements WD/Fa of a cantilever plate with a 
sinusoidally varying point load (P) at the center of the outer 
edge; viscous damping (c/c c ^=0.1) . 


“l 

Plate center 
f=0.5, |=0 

Load point 

Corner point 
a = b = 0*5 

0.0 

0.111 

0.360 

0.329 

0.5 

.148 

.469 

.436 

0.9 

.438 

1.309 

1 . 268 

0.99 

.578 

1.700 

1.669 

1.0 

.576 

1.690 

1.660 

1.01 

.567 

1.662 

1.636 

1.1 

.381 

1.094 

1.093 

1.5 

.093 

.240 

.261 

2.0 

.042 

.088 

.113 

6.1317 

.172 

.404 

.187 

7.8307 

.119 

.129 

.287 

15.6165 

.117 

.200 

.207 


Although the method was demonstrated for homogeneous beams and plates of 
uniform thickness, the energy functionals presented are in more general form, 
capable of accounting for variable cross-sections (i.e., A, I, h) and nonhomo- 
geneous material (E, p) . JWork is currently progressing to extend the method 
further to shell problems. The method will also be developed for analyzing 
the free, damped response of structural elements, and other types of damping 
will also be considered. 
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